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Abstract 

We show that the constituent quark model leads to simple expressions for the 
isoscalar and isovector tensor susceptibilities of the vacuum. The found values are 
negative and of magnitude compatible with QCD-sum-rule parameterizations of 
spectral densities in appropriate L = 1-meson channels. 



Susceptibilities of the vacuum are important quantities of strong interaction physics. They 
directly enter in the determination of hadron properties in the QCD sum rule approach. 
In particular, tensor susceptibilities of the vacuum [1-3] are relevant for the determination 
of tensor charges of the nucleon [4,5]. Widely different values for these susceptibilities, 
all obtained by QCD sum rules techniques, have been reported in the literature [1-3]. In 
this paper we address the issue from a completely different viewpoint, using the concept 
of constituent quarks. It is believed that chiral constituent quark models [6,7] describe 
properly the essential physics in the shallow Euclidean region, with momentum transfers 
of the order of -C A^, where A is the ultraviolet cut off, typically of the order of a 
few hundred MeV. This is strongly supported by the instanton-liquid model [8,9]. In that 
region spontaneous breakdown of chiral symmetry yields a constituent mass for quarks, 
M ~ (250 -r- 400)MeV. Susceptibilities are quantities defined with = 0, hence they seem 
to perfectly fit in the applicability range of the constituent quark model. 

We define the isoscalar and isovector tensor susceptibilities of the vacuum as 
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n,=o = (o|i / d^xr{V'(x)(7'^Xx), v(o)(7^.T^(o)}|o), (i) 

n,=i = (0|z J d^xT{i^{x)a''''T^i;{x),i^{0)a^,T^ij{0)}\0), (2) 

where = (m, d) is the quark field, a^" — |[7^, 7^], and denote Pauh isospin matrices. 
Related definitions are present in the literature: He and Ji [1] introduce n(0) = ^{Uj^q + Ili=i) , 
and Belyaev and Ogancsian [2] consider ni(0) = ^n(O). The main result of this paper is 
that the constituent quark model (CQM) predicts: 

ntT = n?_T = (3) 

9a 9a 



where = 93MeV is the pion decay constant, is the axial charge of the constituent 
quark, and nf^Q ^ denote its isoscalar and isovector anomalous magnetic moments. The 
result (3) is the leading- A^c result, with g'^ and nf^Q i entering as model parameters. The 
expected range for g*^ is 0.75 < g'^ < I, with the lower value naively helping to reproduce 
gA of the nucleon in the non-relativistic quark model. Such lower values are obtained when 
constituent quarks attract in the vcctor-isovector channel [7,10,11]. In the absence of such 
interactions g-^ = 1. It is commonly accepted that the magnitudes of Kj^q Kjj—^ are 
tiny, less than ~ 5%. In constituent quark models these quantities depart from zero when 
tensor interactions among quarks are introduced [10,11]. Using above estimates we find 

- 0.25GeV^ < n^^fi < -0.2GeV^ (4) 



This value is very close to the estimate of Belyaev and Oganesian [2] made in the framework 
of the QCD sum rules: (Hj^i + iUi=o)/2 = -0.2GeV^ 

Now we pass to the derivation of Eq. (3). For our purpose no details of a specific constituent 
quark model are needed, such as the way of introducing the ultraviolet cut off, its value, or 
the value of the constituent quark mass. It is enough that the spontaneous symmetry has 
been broken and quarks are massive. In addition, our calculation is made in the leading-Nc 
approximMion. which is equivalent to the one-quark-loop approximation in the effective 
action. Corrections to this leading result could in principle be included. Wc begin with the 
simplest case of no tensor nor vector interactions. Then Eqs. (1-2) are evaluated by the 
one-quark-loop diagram, which in the momentum space gives 



l-rCQM _ AT- Ar • / 

d^k TrDirac[(|^ + M)a^'^(|^ + M)a^,J _ 
, /• d^k ^ f d^kp. M- 



A 



where A'^c = 3 is the number of colors, Nf = 2 is the number of flavors, and M is the mass of 
the constituent quark. In the last equality we have Wick-rotated to Euclidean momentum 
kE- The integral in Eq. (5) is divergent, as is inherent in effective models, hence a regulator 
is introduced in the integral in Eq. (5), indicated by the subscript A in the integral. There 
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arc many ways of introducing the cut-off, tfie details, fiowever, are irrelevant. The crucial 
point is that on the RHS of Eq. (5) we recognize the well-known result of constituent quark 
models [6,7,10,11] for the square of the pion decay constant: 



Fl = AN, 



Combining Eqs. (5-6) proves Eq. (3) for the case of no vector-isovector {g^ = 1) and 
no tensor (k/=^oi — 0) interaction. Note that our result holds also for the case where M 
depends on the momentum, M — M(k^), as in the case of the instanton- vacuum model 
described in Ref. [6]. 

In the general case interactions in the vector and tensor channels can be present. The 
following terms in the effective chiral SU(2)(8)SU(2) Lagrangian density are relevant for 
our purpose: 



(^(a;)7'^rXa;))^ + (^/'(a;)7'^75TXa;' ' ^ 



x)(7'*Xa;))^ - {i/j{x)a^'''r''i/j{x)y 



(7) 



The signs inside brackets are fixed by the chiral invariance, and the signs outside brackets 
are conventional. The presence of interactions (7) causes rescattering of quarks in vector 
and tensor channels: for instance, the qq pair induced by the tensor current interacts, 
while propagating, by the tensor interactions. The resulting Bethe-Salpeter chain can be 
summed up as a geometric series. Details of this procedure are well known [7] and we give 
no further details here. The effects of finite Gp and Gt are particularly simple for the 
vanishing external momentum, q = 0, which is the case needed here. We find 



^CQM _ 



J(0) 



1 ± 4Gr J(0) 



(8) 



where J(0) = 4iV, MV(A;|- 



M^y. Also, the expressions for the pion decay constant 



and the weak charge of the constituent quark are modified appropriately by the presence 
of vector-isovector interactions [7,10,11], 



1 



1 + 4Gp J(0) • 



(9) 



The anomalous magnetic moments of the constituent quark are equal to [7,10,11] 



7=0,1 



^8GtJ(0)_ 
1 ± 4G'r J(0) 



, or : 1 



7=0,1 



2 1 ± 4G't J(0) ' 



(10) 



With the help of relations (9,10) we rewrite Eq. (8) as Eq. (3), completing the proof. 
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In the remaining part of this paper we argue that the sign and size of n/=o and n/=i, 
as given by Eq. (4), agree with expectations based on dispersion relations and simple 
models of spectra in appropriate channels. Let us begin with Ili=i, for which the rele- 
vant intermediate states have quantum numbers of p {I^{J^^) = 1^(1 )) and 6i(1235) 
(/'^(J^'") = 1^(1^^)). With the help of the spectral decomposition we can write the imag- 
inary part of the correlator n/=i(g) = i{0\ J d'^x e^'^'^T{ip{x)a^''T^ip{x),tp{0)af^uT^tjj{0)}\0) 
as 



Imn,=i(g) = TT 5: 6{q' - 5:(0|^(0)a^V=^V^(0)|p; g, A) (p; g, \\m<y,vT^m^) + 

P A 

+ TT 5: (5(g2 - ml) 5:(0|^(0)a'^V^(0) |6; g, A) (6; g, A|^(0)a^,r3^(0)|0), (11) 



where X]p,fe denote the sum/integral over all physical states with quantum numbers of 
neutral p and 6i, quantities rup and denote the invariant mass of these states, and A are 
their polarizations. The Lorentz structure of the matrix elements consistent with quantum 
numbers of the state isQ 

(0|^(0)a'^V3^(0)|p; g, A) = i^fM^) {e^.q" - e^q^) , (12) 
(0|V^(0)cr^V3^(0)|6; g, A) = iV2h{q^)e^-^P{e^)^qp, (13) 

where is the polarization vector. Using the identity X^a ^aI^a)* — ~9'^^ \q^(f jc^ we find 
that 

$:(0|V^(0)a'^V^(0)|p; g, A)(p; g, A|V^(0)a^,r3^(0)|0) = -Vlq^\f ,{q^)\\ (14) 

A 

$:(0|V^(0)a'^VV^(0)|6; g, A)(6; g, A|^(0)a^,r3V^(0)|0) = ^\2q^\M)\\ (15) 

A 

and, consequently, Eq. (11) can be written as 

Imn,=i(g) = -127r^5(g2 - mj)mj|/,(mj)p + 127r^5(g2 - ml)ml\U{ml)\\ (16) 



Note the negative sign of the p contribution, and the positive sign of the h\ contribution. 
The dispersion relation gives 

nj=i(g) = — / ds ' , 17 

IT J s — q'^ — %o 



where in general the spectral densities collect the contribution from all poles (and cuts) in 
the appropriate channel: cr(s) = Yja^{s—m\)m'^^\fa{rn?^)\'^. Dispersion relation (17) requires 
no subtraction in the chiral limit. This follows from the fact that in QCD the asymptotic 

^ The factors of \/2 in Eqs. (12-13) make our expressions consistent with the definition 
(0|u(0)a'^-d(0)|p+; (7, A) = ifp{q^) [e^^q^ - e'^q'^) of Ref. [12]. 
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forms at large Euclidean momenta (g^ — oo) of tensor correlators in the p and bi channels 
are equal [13]. Indeed, the tensor correlator 



can be decomposed as as follows in structures involving the contributions of p and bi states: 

U^jfiq) = -4T(q"%(g) + 4T(7f ^n,(g), (19) 
where the factors of 4 are conventional^. The tensor structures are defined as 

= - q^q^g'^'' + q^q'^g''^ - q'^q'^g^^) , (20) 

Zq 

Ti;^f = -Tl'Y + 1 {a^'a'"' - g'^-'g^") , (21) 

and satisfy Tj^^^^^^ = Tj^^^^^^ = 3. In the deep Euclidean region QCD gives [12,13][3- 



nr fe) = -If j (l ± l°g ^ + 0(a.)\ + (22) 



where m is the current quark mass, which is neglected. The form (22) implies that at large 
s we have asymptotically (Jp{s) ~ (Jb{s) ~ s/(87r), and the difference crp(s) — crb{s) ~ 1/s. 
This means that the unsubtracted dispersion relation (17) holds, and we obtain immediately 



oo 



n,_, = i^/d.^^^i^l^^. (23) 

71 J S 


Repeating the above analysis for the isovector channel, where the relevant physical states 
have the quantum numbers of (/^(J^^) = 0-(l — )) and /ii(1170) (/^(J^^) =0^(1+-)), 
with conventionally the same factors as in Eqs. (12,13,19,22) gives 



oo 



n,-o = lg/d /^(^)-"-^^) . (24) 
n J s 





Note that in Eqs. (23-24) there is cancelation between spectral densities corresponding to 
opposite parity states. 



^ We note that the constituent quark model gives ll^^^(0) + 11^^^ (0) = 0, in agreement with 
the general requirement due to the absence of zero-mass states in the tensor channel [12,14]. 
^ Ref. [13] has an error in the overall factor in Eq. (A3). Their tensors in Eqs. (2.3) relate to ours as 
follows: = ir(_), ^ -lr(+), and Eq. (A3) should read: n(g) = 24P(i)n- + 24P(i)n+, 
in order to agree with our Eq. (19) 
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Now, following Refs. [1,2,12], we can estimate the contributions to Eq. (23-24). We expect 
that the following pole + continuum parameterization of the spectral densities should give 
at least a correct order-of-magnitude: 



aa{s) = nmlfJis - ml) + ^e{sa - s), (25) 



where a labels the channel. Substituting Eq. (25) into Eqs. (23-24) gives 



n,=o = + -^(5, - s^)), (26) 

n,=i = -12(f^^-f,2 + ^(.,-.,)) (27) 

Let us begin with the last term in Eq. (27). The Sf, parameter has been obtained from QCD 
sum rules in Refs. [2,12,13] giving ~ (2.3 -j- 3)GeV^. Using derivative coupling Ref. [15] 
reports ~ 2.5GeV^. The QCD sum rule fit for the p meson spectrum in the tensor chan- 
nel gives [12,14,16] Sp ~ (1.3 ^ 1.7)GeV^. Thus we have Sf, - Sp - 0.6 ^ 1.7GeV^ and the 
continuum contribution to I\.i=i is of the order of — O.lGeV^ to — 0.25GeV^. The hi pole 
contribution is, according to Ref. [2,12], 12f^ ~ 0.4GcV^, and the p pole contribution is, 
according to Ref. [12], 121^ ~ 0.3GeV^ according to [16], 121^ ~ 0.35GeV^ and according 
to [14], 12f2 - 0.5GeV^ Taking a "global" average yields n/=i - (-0.1 ^ -0.25)GeV^ 
Note that the order of magnitude of all contributions to Eq. (27) is the same as in Eq. (4), 
indicating our constituent-quark-model value has the expected size. The sign of Eq. (4) 
shows that the p channel "wins" over the hi channel in Eq. (27). In fact, as already pointed 
out in Ref. [1], cancelations in Eqs. (26-27) make any precise estimates of tensor suscepti- 
bilities based on models of spectra difficult if not impossible. Note, however, as pointed out 
by Belyaev and Oganesian [2], that Ref. [1] does not include the continuum contribution 
in Eq. (27), which causes their value to be much smaller, and of opposite sign compared 
to Ref. [2] and to our estimate (4). The value for 11(0) obtained by Kisslinger [3] from the 
QCD sum rules for three-point functions are of opposite sign and large in magnitude. In 
view of these controversies, direct model estimates, such as that of Eq. (4) gain significance. 

We point out that the mentioned QCD sum rule calculations do not distinguish between 
the isoscalar and isovector channels. Surprisingly, it is supported by or result (3), since, 
due to the smallness of the quark anomalous magnetic moments, IIj^!^ is practically equal 
to IIj^L^- A corollary of this observation is the approximate relation 
/o°° ds{(Th - (7^)/s ~ /o°° ds{ai - (Tp)/s. 

We conclude with several comments: 

(1) The obtained values for tensor susceptibihties of the vacuum have direct relevance for 
the calculations of tensor charges of the nucleon [1,4,5]. 

(2) The simple result (3) and the following numerical estimate (4) are leading-A'^c re- 
sults. Formally suppressed effects, such as meson (especially pion) exchange, may 
be important, as is the case e.g. in the calculation of the quark condensate in the 
Nambu-Jona-Lasinio model [17]. 
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(3) Quark-model relations, such as our Eq. (3), bear some similarity to the relations 
obtained by Leutwyler [18] and Mallik [19] on the basis of the SU(6) symmetry of 
meson wave functions. In fact, more relations of the type of Eq. (3) can be obtained 
in our approach [20]. 

(4) Our relations of susceptibilities to the L = 1 meson spectra, Eqs. (23-24), are obtained 
via dispersion relations. This is in the sprit of QCD sum rules, where one compares the 
deep Euclidean region to the physical region via dispersion relations. In our method 
we compare the shallow Euclidean region to the physical region, following the idea 
described in Ref. [21]. That way we can infer useful information concerning meson 
excitations without leaving the validity range of the constituent chiral quark model. 

We thank Rusko Ruskov for many discussions on vacuum susceptibilities in QCD sum 
rules. 
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